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Abstract 

We have derived long series expansions of the percolation probability for 
site, bond and site-bond percolation on the directed triangular lattice. For the 
bond problem we have extended the series from order 12 to 51 and for the 
site problem from order 12 to 35. For the site-bond problem, which has not 
been studied before, we have derived the series to order 32. Our estimates of 
the critical exponent (3 are in full agreement with results for similar problems 
on the square lattice, confirming expectations of universality. For the critical 
probability and exponent we find in the site case: qc = 0.4043528 it 0.0000010 
and (3 = 0.27645 ± 0.00010; in the bond case: qc = 0.52198 ± 0.00001 and 
13 = 0.2769 ± 0.0010; and in the site-bond case: qc = 0.264173 ± 0.000003 and 
/3 = 0.2766 lb 0.0003. In addition we have obtained accurate estimates for the 
critical amplitudes. In all cases we find that the leading correction to scaling 
term is analytic, i.e., the confluent exponent A = 1. 



PACS numbers: 05.50.+q, 02.50.-r, 05.70.Ln 



1 Introduction 



In an earlier paper (Jensen and Guttmann 1995) we reported on the derivation and 
analysis of long series for the percolation probability of site and bond percolation on 
the directed square and hexagonal lattices. In this paper we extend this work to site, 
bond and site-bond percolation on the directed triangular lattice. We refer to our 
earlier paper for a more general introduction to directed percolation and its role in the 
modelling of physical systems. In directed site percolation each site is either present 
(with probability p) or absent (with probability q = I — p) independent of all other 
sites on the lattice. Similarly for bond percolation each bond is absent or present 
independently of other bonds. Finally in site-bond percolation both sites and bonds 
may be absent or present with equal probability, but again with no dependency on 
any other sites or bonds. Two sites in the various models are connected if one can 
find a path, respecting the directions indicated in Figure through occupied sites, 
bonds or sites and bonds, respectively, from one to the other. When p is smaller than a 
critical value Pc all clusters of connected sites remain finite, while for p > Pc there is an 
infinite cluster spanning the lattice in the preferred direction. The order parameter of 
the system is the percolation probability P{p) that a given site belongs to the infinite 
cluster. This quantitiy is strictly zero when p < pc and changes continuously at pc- 
For p > Pc the behaviour of P{p) in the vicinity of Pc may be described by a critical 
exponent P, 

P{p) (x{p-pcf, p-*pt- (1) 

The bond problem was originally studied by Blease (1977) who calculated a series 
to 12th order. For the site problem De'Bell and Essam (1983) derived the series 
to 12th order. The site-bond problem has, at least to our knowledge, never been 
studied before. Our main motivation for doing so in this paper is to obtain further 
independent estimates of the critical exponent (3. Using the finite-lattice method and 
the extrapolation technique of Baxter and Guttmann (1988) we have extended the 
series for the bond problem to order 51, for the site problem to order 35 and derived 
the series for the site-bond problem to order 32. The site and bond problems have also 
been studied by Essam et al. (1986,1988), who derived series expansions for moments 
of the pair connectedness. 

2 The finite- lattice method 

We wish to derive a series expansion for the percolation probability on the directed 
triangular lattice oriented as in Figure ||. In this figure we have numbered the various 
levels or rows of the lattice according to which sites can be reached by a path of 
minimum length A^ — 1 starting at the origin O. In other words all sites in the iVth row 
can be reached in — 1 steps but not in — 2 steps. Note that a path going through 
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Figure 1: The directed triangular lattice with orientation given by the arrows. The rows 
are labelled according to the text. 

a given site can only reach the part of the lattice shown in Figure |1| below the origin 
O. This suggests that one should look at the following finite-lattice approximation to 
P{q), namely the probability Pn{.(1) that the origin is connected to at least one site 
in the iVth row. Since we are in the high density region we have chosen to use the 
expansion parameter q rather than p. Pn{(i) is a polynomial with integer coefficients 
and a maximal order determined by the total number of sites and/or bonds on the 
finite lattice. 

By the method used for the square lattice problems (Bousquet-Melou 1995) one can 
prove, mutatis mutandis, that the polynomials -P/v(g) converge to P{q). Indeed we 
may consider P{q) = limTv^oo -P/v(?) to be a more precise definition of the percolation 
probability. More importantly, however, from a series expansion point of view, for 
the site and site-bond problems the first + 1 coefficients of the polynomials -P/v(q') 
are identical to those of P{q). In the case of bond percolation the agreement extends 
through the first 2N + 1 coefficients. 

2.1 Specification of the models 

To calculate the finite-lattice percolation probability Pniq) "we associate a state aj 
with each site, such that o"-,- = 1 if site j is connected to the A^th row and aj = — 1 
otherwise. We shall often write +/— for simplicity. Let /, c and r denote the sites 
connected to a site t from the row above, as in Figure |l[ We then define the weight 
function W{at\cri, o^, Or) as the probability that the top site t is in state Of, given that 
the lower sites /, c and r are in states a;, cXc and 0"^, respectively. As for the square 
lattice (Bidaux and Forgacs 1984, Baxter and Guttmann 1988) we then have 
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{a} t 



(2) 



where the product is over all sites t of the lattice above the iVth row. The sum is over 
all values ±1 of each at, other than the topmost spin ai which always takes the value 
+1. The spins in the iVth row are fixed at +1, and Pn{<i) is calculated as the sum 
over all possible configurations of the probability of each individual configuration. 

The weight functions W are calculated as follows. Obviously, W{—\ai,ac,(Jr) = 
1 — W{+\ai,ac,crr)- The remaining weights are easily calculated by considering the 
possible arrangements of states and sites and/or bonds. W{+\—,—,—) = because 
the top site is connected to the A^th row if and only if at least one of its neighbours 
below is connected to the A^th row. All the remaining weights for the site problem 
equal I — q because the top-site has to be occupied in order to be connected to the 
A^th row. Let us next look at the remaining bond weights. ^^'^(+1+, +,+) = 1 — 
because the only bond configuration not allowed is all three bonds absent, which has 
probability W^{+\+, +, -) = W^{+\+, -, +) = ^"^(+1-, +, +) = 1 - because 
the bond to the — state can be either present or absent (probability 1) while among 
the remaining bonds only the configuration with both bonds absent (probability g^) 
is forbidden. Finally, W^{+\+,-,-) = W^{+\-,+,-) = W^{+\-,-,+) = l-q 
because the bond to the + state has to be present, which happens with probability 
p = 1 — q, while the other bonds can be either present or absent. For the site-bond 
problem we find that W^^{+\ai,ac,<Jr) = (1 — q)W^ {+\ai, ac, <Jr) because if the top 
state is +1 the top site has to be present. 

2.2 Series expansion algorithm 

Computer algorithms for the calculation of -P/v(q') are readily found. These are basi- 
cally implementations of the transfer matrix technique. The general features of these 
algorithms were described in our earlier paper (Jensen and Guttmann 1995), to which 
we refer for further details. The sum over configurations is performed by moving a 
boundary line through the lattice. For each configuration along the boundary line one 
maintains a (truncated) polynomial which equals the sum of the product of weights 
over all possible states on the side of the boundary already traversed. The boundary 
is moved through the lattice one site at a time. The calculation of -P/v(g) by this 
method is limited by memory, since one needs storage for 2^ boundary configurations. 
However, as was the case with the square lattice, this problem can be circumvented 
by introducing a cut into the lattice. For each fixed configuration of states on this 
cut one evaluates the lattice sum P^{q) and gets Pwiq) = J2cPNil) ^^e sum over 
all configurations of the cut. By placing the cut appropriately, the growth in memory 
requirements can be reduced to 2^^^. 

In Figure ^ we show the triangular lattice with a cut marked by filled circles. In the 
algorithm the cut is used as a pivot line by the boundary line which traverse the lattice. 
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Figure 2: The directed triangular lattice with orientation given by the arrows. The sites 
with fixed states along the pivot line are marked by filled circles. The open circles mark 
one particular position of the boundary line during the traversing of the lattice. 

We start by building up the first row at the base CL of the lattice. We then build 
up the part of lattice above the cut from row CL to row EL'. Next the boundary line 
expands along the line-piece ES until it reach the position ESL" and the last site (at 
L") is flipped to the other site of the top-most triangle (after this the boundary line 
is in the position marked by the open circles). Then we work our way down the right 
side of the lattice past R to position ESB. Finally the boundary line is moved down 
along the line-piece SEC after which the whole lattice has been build up. This process 
is then repeated for each configuration of the cut. Since the calculations for different 
cut-configurations are independent of each other this algorithm is perfectly suited to 
take full advantage of massively parallel computers. 

Using this algorithm we calculated Piv(?) for N < 23 for the bond and site-bond 
problems. The integer coefficients of -PAr(g) become very large so the calculation was 
performed using modular arithmetic (see, for example, Knuth 1969). Each run with 
N = 23, using a different moduli, took approximately 70 hours for the bond problem 
and 55 hours for the site-bond using 50 nodes on an Intel Paragon. For the site problem 
the weights only depend on whether or not there are any H-'s among the neighbours 
of the top-most site. As was the case for the square site problem this may be used 
to sum over many configurations of the cut simultaneously (see Jensen and Guttmann 
(1995) for further details). This allowed us to calculate PN{q) for N < 25. Each run 
for N — 25 took about 85 hours using 50 nodes. 
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3 Extrapolation of the series 



As mentioned, the coefficients of tlie polynomials Pn{(i) = J2m>o 0'N,m(f^ will generally 
agree with those of the series for P{q) = J2m>o o-mff^ up to some order, A^, determined 
by A^, but depending on the specific problem. In the case of directed bond percolation 
on the square lattice Baxter and Guttmann (1988) found that the series for P{q) could 
be extended significantly by determining correction terms to Pn{q)- Let us look at 

PN-PN+l = q^J2^''dN,r (3) 
r>0 

then we call (ijv,r = (ijyN+r ~ ^N+iN+r correction term. If formulas can be 

found for d^.r for all r < A' then, using the series coefficients of PnIq), one can extend 
the series for P{q) to order N + K since 

r 

^N+r ~ ^N,N+r ~ X] ^N+r-m,m (4) 
m=l 

for all r < K. That this method can be very efficient was demonstrated by Baxter 
and Guttmann, who identified the first twelve correction terms for the square bond 
problem, and used P29{q) to extend the scries for P{q) to order 41. To really appreciate 
this advance one should bear in mind that the time it takes to calculate PN{q) grows 
exponentially with A", so a direct calculation correct to the same order would have 
taken years rather than days. In the following we will give details of the correction 
terms for the various directed percolation problems on the triangular lattice. 

3.1 The site problem 

For the site problem the coefficients of PN{q) agree with those of P{q) to order A". In 
this case the first correction term is very simple as d^^o = 2 for A" > 2, i.e., the first 
correction term is simply a constant. For the second correction term o?7v,i we find the 
following sequence 

0,0,3,18,32,50,72,98,... 
It is thus immediately clear that 

dN,i^2N^ for A^>3. (5) 

Note that for convenience we assume that the sequence starts from A^ = 0. And indeed 
we find that for A^ > r + 1, dN,r can be expressed as a polynomial in A^ of order 2r. Wc 
have been able to calculate these polynomials for the the first 10 correction terms. It 
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k/r 1 2 3 4 5 6 7 8 9 
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-345600 


-65318400 
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-6720 


662400 


86728320 
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160 


-2256 


-299136 


-54616320 


3 




-12 


-456 


-5592 


-155040 


29156640 


4 
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112 


6968 


262400 


3721088 


5 






-72 


-4680 


-211440 


-13781520 


6 






16 


1016 


117072 


9766720 


7 








-288 


-35760 


-3900960 


8 








32 


6000 


1183584 


9 










-960 


-222000 


10 










64 


28000 


11 












-2880 


12 












128 



13 
14 
15 
16 
17 
18 



-15850598400 


-2984789606400 


-539895767040000 


15417077760 


3039204188160 


681914690150400 


-10042993152 


-2801552624640 


-758646639912960 


6930400512 


1683396497664 


492391103938560 


-1895857152 


-641242189440 


-236796916234752 


-275292864 


183056948928 


80349078951936 


775939360 


32888441824 


-13942053553664 


-484442784 


-52810790592 


-3002221192320 


180360704 


27746932192 


4062978111936 


-46002432 


-9468263616 


-1860005271168 


8946336 


2268003232 


567526218432 


-1175328 


-405615168 


-128527251840 


112448 


55739936 


21947992384 


-8064 


-5494272 


-2918143872 


256 


406784 


301743168 




-21504 


-23270400 




512 


1362432 






-55296 






1024 



Table I: The coefficients cjr in the extrapolation formulas Eq. (H) for the site problem. 

turns out that it is useful to pull out a factor l/(r!(r + 1)!) and express the correction 
terms as 

1 2r 

r!(r + 1)! 

This ensures that the coefficients in the extrapolation formulas are integers. We 
have listed these coefficients in Table H 

Obviously since these formulas are correct for N > r + 1 and we have calculated Pn{(i) 
for < 25 we did not have enough terms in the correction sequences to calculate all 
the coefficients in these polynomials for the largest values of r. However, from the 
table of coefficients, it is immediately clear that cfj' = . And in general we found 
that (?^~'^ 12^^^ is a polynomial in r of order 2m 



or-l-l 2m 

C^-'" = — ^ E ^m^', (7) 



-4)m^l 



fe=0 



where the prefactor has been chosen so as to make the leading coefficients particularly 
simple. In Table |T| we have listed the coefficients 6^ for the first six polynomials. 

This time we note that 6^ = 3"^. And indeed as before we find that 6^~''/3™' is a 
polynomial in m of order 2j. In particular we have. 
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11 












12474 
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729 



Table II: The coefficients 6,^ in the extrapolation formulas Eq. ([7|) for the site problem. 

^2rn-i _ 3^^(17/27+ 10/27m), 

and 

^2m-2 ^ 3'»m(ioi5/486 - 5137/1458m + 332/243m2 + 50/72W). 

So when calculating the extrapolation formulas Eq. we first used the sequences 
for the correction terms to predict as many polynomials as possible. When we ran 
out of terms we then predicted as many of the leading coefficients from Eq. (0) as 
possible. This in turn allowed us to find more extrapolation formulas, which we could 
use (together with the formulas for &^~-' ) to find more of the formulas for c^^"™". And 
so on until the process stopped with the 10 extrapolation formulas we listed above. 

Using the 10 extrapolation formulas and P2h{<l) we extended the series for P{q) through 
order 35. The resulting series is listed in Table |in[ 



3.2 The site-bond problem 

For the site-bond problem the coefficients of Pn{(i) agree with those of P{q) to order 
A^. In this case the correction terms are very similar to those of the site problem. In 
particular we find that d^^ = 12 and in general d^^r is a polynomial in of order 2r, 

or 2r 
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-111307 
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-255236 
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-590543 
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-1 
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-1362919 
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-2 
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-5 
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-182 
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-539266908 
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-393 
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-1200140540 


12 


-853 


30 


-3087966932 


13 


-1887 


31 


-6454135923 


14 


-4208 


32 


-18281313306 


15 


-9445 
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-33072764132 


16 


-21350 


34 


-114854030873 


17 


-48612 


35 


-145978838818 



Table III: The coefficients a„ in the series expansion of the percolation probability P{q) 
for directed site percolation on the triangular lattice. 

We have identified the first 9 correction terms for the site-bond problem and have 



listed the coefficients in the extrapolation formulas in Table IV. 



From this table it is immediately clear that the coefficient of the leading order c^^ = 
3 X 4^. As in the site case we find that c^^~^ / A^'^^ is a polynomial in r of order 2m. 



4r+l 



2m 



2r—m \ ' 7fc „k 



E btr\ (9) 



where the prefactor has been chosen so as to make the leading coefficients particularly 
simple. In Table we have listed the coefficients for the first six polynomials. 

In this case 6^ = 3"^+^ and b'^"^ = 3'"'^^m(10/27 — 16/27m), which, using the same 
procedure as before allowed us to find the first 9 extrapolation formulas. From P23{q) 
we were thus able to extend the series for P{q) through order 32. The resulting series 



is listed in Table VI 



3.3 The bond problem 

For the bond problem the coefficients of Pn{(i) agree with those of P{q) to order 2N. 
In this case the first correction term is more complicated. For the first correction term 
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Table IV: The coefficients in the extrapolation formulas Eq. (|) for the site-bond 
problem. 
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Table V: The coefficients 6^ in the extrapolation formulas Eq. (^) for the site-bond 
problem. 
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Table VI: The coefficients a„ in the series expansion of the percolation probability P{q) 
for directed site-bond percolation on the triangular lattice. 

d^fl we find the following sequence 

1,3,9,27,83,263,857,... 

which we have identified as 

rfiv,o = '2Cn - 1. (10) 

where Cn = {2N)\/{{N + 1)!A^!) are the Catalan numbers, which also occur in the 
correction terms for the square bond problem. In general we find that for r < 4 the 
correction terms are given, for > r — 2, by the formulas 

r+l / N \ 1 2r 

dN,r = E^rCN^'^-l + Eb'r( , ] CN + -T;^ ^ ^'N' ■ (H) 
k=l k=l \ / ' ■' ■ k=0 

We have listed the coefficients a^, 6^ and of these extrapolation formulas in Table [VI 1| . 
We note that as in the previous problems the leading coefficients are quite simple, 

a;+i = (-l)'-2a+i, K = 2, and c^/ = -C,.. 

These 5 extrapolation formulas and P2'i{<l) allowed us to extend the series for P{q) 
through order 51. The resulting series is listed in Table |V111[ 
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Table VII: The coefficients a^, 6^ and in the extrapolation formulas Eq. (|TI|) for the 
bond problem. 
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Table VIII: The coefficients a„ in the series expansion of the percolation probability P{q) 
for directed bond percolation on the triangular lattice. 
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4 Analysis of the series 



We expect that the series for the percolation probabihty behaves hke 

P{q) ~ A{1 - q/q,)^[l + a^{l - q/q,)^ + ...], (12) 

where A is the critical amplitude, A the leading confluent exponent and the . . . repre- 
sents higher order correction terms. In the following sections we present the results of 
our analysis of the series which include accurate estimates for the critical parameters 
gc, A and A. For the most part the best results are obtained using Dlog Fade (or in 
some cases just ordinary Fade) approximants. A comprehensive review of these and 
other techniques for series analysis may be found in Guttmann (1989). 



4.1 qc and jS 

In Table [TX] we list various Dlog Fade approximants to the percolation probability series 



for directed site percolation on the triangular lattice. The defective approximants, 
those for which there is a spurious singularity on the positive real axis closer to the 
origin than the physical critical point, are marked with an asterisk. Most higher-order 
approximants yield estimates around the values qc = 0.4043528 and f3 = 0.27645, with 
very little spread among the approximants. Opting for a conservative error estimate, 
it seems appropriate to estimate that the critical parameters lie in the ranges, qc = 
0.4043528(10) and (3 = 0.27645(10), where the figures in parenthesis indicate the 
estimated error on the last digits. 

The results of the analysis of the series for the bond problem are listed in Table ^. 
In this case the spread among the various approximants is quite substantial, there 
appears to be a marked downward drift in the estimates for both qc and /?, and the 
estimates do not settle down to definite values. It does however seem likely that the 
true critical parameters lie within the estimates: qc = 0.52198(1) and (3 = 0.2769(10). 

The analysis of the series for the site-bond problem yields the results in Table 0. 
Again we see a downward drift in the estimates for both qc and (3 though the estimates 
are somewhat more stable than in the previous case. We estimate that the true critical 
parameters lie within the ranges: qc = 0.264173(3) and P = 0.2766(3) 



4.2 The critical amplitudes 

We can estimate the critical amplitude A by evaluating Fade approximants to G{q) = 
(Qc — q)P^^^^ at qci since it follows from the leading critical behaviour in Eq. ([T2|) that 
G{qc) ~ A'^/^qc- This prodecure works well but requires knowledge of both and 
(3. As we have just shown, we know both and (3 very accurately for the triangular 
site series. We estimated A using values of qc between 0.4043524 and 0.4043534 and 
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Tabic IX: Dlog Pade approximants to the percolation series for directed site percolation 
on the triangular lattice. 
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Table X: Dlog Pade approximants to the percolation series for directed bond percolation 
on the triangular lattice. 
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Table XI: Dlog Pade approximants to the percolation series for directed site-bond perco- 
lation on the triangular lattice. 

values of /5 ranging from 0.2764 to 0.2765. For each {Qc, pair we calculate A as the 
average over all [N + K, N] Pade approximants with K = 0, ±1 and 2N + K > 25. 
The spread among the approximants is minimal for Qc — 0.4043527, P — 0.27645 
where A — 1.581883(5). Allowing for values of Qc and P within the full range we get 
A = 1.5819(4). 

For the bond problem we used values of Qc from 0.52196 to 0.52121 and P from 0.2763 
to 0.2773 averaging over Pade approximants with 2N+K > 40. In this case the spread 
is minimal for = 0.521985, P = 0.2767 where A = 1.48584(2). Again allowing for a 
wider choice of critical parameters we estimate that A — 1.486(6). 

For the site-bond series we restricted Qc to lie between 0.264170 and 0.264176 and P 
between 0.2763 to 0.2768 using all approximants with 2A'" + K > 25. The minimal 
spread occurs at Qc = 0.264173, P = 0.2766 where A = 1.477393(4). A wider choice 
for Qc and P leads to the estimate A — 1.477(1). 

4.3 The confluent exponent 

We studied the series using two different methods in order to estimate the value of 
the confluent exponent. In the flrst method, due to Baker and Hunter (1973), one 
transforms the function P, 

n oo 

P{q) = Y: A{1 - q/qc)-'' = E ^n^'^ (13) 

i=l n=0 
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into an auxiliary function with simple poles at 1/Aj. We first make the change of 
variable q = qdl — e~'^) and find, after multiplying the coefficient of by k\, the 
auxiliary function 

N oo N' A 

(14) 



i=l k=0 i=l '^'^^ 



which has poles at C = with residue — Aj/Aj. The great advantage of this method 
is that one obtains simultaneous estimates for many critical parameters, namely, /3 (the 
dominant singularity), A (the sub-dominant singularity), and the critical amplitudes 
(the residues at the singularities), while there is only one parameter Qc in the trans- 
formation. Unfortunately this method does not appear to work well for this problem. 
For the site problem we find that the transformed series generally yields poor esti- 
mates for f3 and no estimates for the confluent exponent. For the bond and site-bond 
problem the situation is somewhat better. In Table |X11| we have listed estimates for 



the critical parameters obtained from various Fade approximants to the Baker-Hunter 
transformed series, using the values Qc = 0.52198 for the bond series and qc = 0.264173 
for the site-bond series. 
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13 


13 


0.27342 
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1.01091 


2.16997 


14 


15 


0.27661 


1.47745 


0.99950 


2.08954 


15 


15 


0.27828 


1.48182 


0.96056 


1.91013 


15 


16 


0.27659 


1.47728 


1.00194 


2.10641 



Table XII: The critical exponent /3, confluent exponent A and critical amplitudes A and 
qa obtained from [N,M] Pade approximants to the Baker-Hunter transformed series for 
the bond and site-bond problems. 
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It should be noted that, obviously, all approximants yield estimates for the critical 
parameters. However, we have discarded many approximants from the table because 
we beheve the results to be spurious. For all the discarded approximants we found 
that the amplitude of the confluent term was of order zero and generally the estimate 
for f3 was very far from the expected value. Among the remaining approximants we 
clearly see that the favoured value of the confluent exponent is A = 1. We also note 
that the amplitude estimates are in full agreement with those of the previous section. 

In the second method, due to Adler et al. (1981), one studies Dlog Fade approximants 
to the function where 



F{q)=l3P{q) + {q,-q)dP{q)/dq. 

The logarithmic derivative of F{q) has a pole at qc with residue /? + A. We evaluate 
the Dlog Fade approximants for a range of values of qc and (3. In Table pClli| we have 



listed the estimates for A obtained by averaging over all [A^, N + K] approximants for 
a few values of (3 with fixed at the central value of our estimate range. For the 
site and site-bond problem we used all approximants with 2N + A^ > 25 and for the 
bond problem all approximants with 2N + K > 40. This analysis clearly indicates 
that A ~ 1 and thus that there is no sign of any non-analytic corrections to scaling. 
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0.98587 


0.27630 


0.97076 
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1.03471 
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Table XIII: Estimates for the confluent exponent A from the transformation due to Adler 
et al. (1981) for various values of [3 at the critical point qc. 
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T bond 


0.52198(1) 


0.2769(10) 


1.486(6) 


0.521971(5) 


1.4841(2) 


45 


T site 


0.4043528(10) 


0.27645(10) 


1.5819(5) 


0.4043523(3) 


1.58183(2) 


30 


T site-bond 


0.264173(3) 


0.2766(3) 


1.477(1) 


0.264170(4) 


1.4765(3) 


25 


S bond 


0.3552994(10) 


0.27643(10) 


1.3292(5) 
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1.32920(1) 


45 


S site 


0.294515(5) 


0.2763(3) 


1.425(1) 


0.294518(3) 


1.42588(4) 


30 


H bond 


0.177143(2) 


0.2763(2) 


1.106(1) 


0.177144(2) 


1.1064(3) 


30 


H site 


0.160067(5) 


0.2763(4) 


1.167(1) 


0.160069(2) 


1.1680(3) 


30 



Table XIV: Estimates of critical parameters for the three problems on the triangular (T) 
lattice studied in this paper and for the site and bond problems on the directed square (S) 
and honeycomb (H) lattices. See the text for explanation of the biased estimates. 



5 Conclusion 



In this paper we have presented extended series for the percolation probability for 
site, bond and site-bond percolation on the directed triangular lattice. The analysis 
of the series leads to improved estimates for the percolation threshold and the order 
parameter exponent (3. In Table XIV we summarise the critical parameter estimates 
for the percolation probability for the three problems on the triangular lattice studied 
here and the problems studied in our earlier paper. The estimates for Qc = 1 — Pc for 
the triangular bond and site problems are in excellent agreement with those obtained 
by Essam et al. (1986, 1988), = 0.40437(7) and = 0.521975(7), respectively. The 
estimates for (3 clearly show, as one would expect, that all the models studied in this 
and our earlier paper belong to the same universality class. The unbiased estimates 
for f3, derived in the manner described in the previous section, for the triangular 
site and square bond cases are in excellent agreement and have small error bars (we 
emphasize once more that our error estimates are conservative). This leads us to 
belive that an improved estimate f3 = 0.27644(3) is reasonable. We used this highly 
accurate estimate to obtain the biased estimates in Table pClV| as follows. First we 
formed the series for P{q)~^/^ using (3 = 0.27644. This series has a simple pole 
at qc which can be estimated from ordinary Fade approximants. By averaging over 
all [A^, N + K] approximants with K = 0, ±1 and 2N + K > Nmin we obtained 
the biased estimates for qc the error-bars are basicly twice the spread among the 
approximants. We then used the biased estimate for q^ (with [3 as before) to obtain 
the biased estimates for the amplitudes using the procedure decribed in the previous 
section. As previously noted (Jensen and Guttmann 1995), there is no simple rational 
fraction whose decimal expansion agrees with our estimate of (3. Given that this 
model is not conformally invariant, and that the expectation of exponent rationality 
is a consequence of conformal invariance, it is perhaps naive to expect otherwise. It is 
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nevertheless true that there is a widely held - if imprecisely expressed - view that two 
dimensional systems should have rational exponents. More precise numerical work such 
as the recent estimation of the longitudinal size exponent z/|| (Conway and Guttmann 
1994) of directed animals and the present calculation, supports the conclusion that 
the critical exponents for these models should not be expected to be simple rational 
fractions. Finally note that none of the series show any evidence of non-analytic 
confluent correction terms. This provides a hint that the models might be exactly 
solvable. 
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Appendix: The first extrapolation formulas 



In this appendix we shall calculate the first correction term(s) dM,r for the various 
problems we have studied in this paper. In the following we rely heavily on the work 
of Bousquet-Melou (1995) and we shall represent the directed percolation models in 
terms of directed animals. By a directed site (bond) animal A we simply understand 
any finite set of connected sites (bonds) starting at the origin O in Figure The 
area (or size) \A\ of an animal is the number of sites in the animal and the perimeter 
p{A) is the number of unoccupied sites (bonds) with a nearest neighbour in A. The 
height h of an animal is the last row to which the animal extends, i.e., there is at least 
one occupied site in row h belonging to A but none in row h + 1. The percolation 
probability, for the site and site-bond cases, is 

where A denotes the set of animals on the lattice. For bond percolation the power 
of (1 — q) in the above equation is \A\. The difference stems from the assumption 
that for site percolation the origin is occupied with probability 1. In analogy with the 
finite-lattice formulation we define subsets of ^ as the set of animals of height less 
than N . It follows that 

PN{q) = l- E /^^^(l-g)l^l-\ (A.2) 

and 

PM{q) - Pn+M = I - E /^^Hl-g)l^l-^ (A.3) 

AeAN+i\AN 

It should be noted that in the site and site-bond cases the polynomials PN^q) defined 
above are identical to the polynomials Pn+i{(i) from Section 2. From Eq. ( |A.3|) it is 
immediately clear that P/v and Pn+i agree up to an order determined by the animals 
in ^7v+i\^Ar with the smallest perimeter. In our cases is simply proportional to 
N and the polynomials PN{q) therefore have a formal limit Poo{q) which we identify 
as the percolation probability P{q). By expanding Eq. ( |A.3| ) one gets a very useful 
expression for the correction terms 



r>0 r>0 k=OA&ANk V ^ / 



where AN,k = {A e ^7v+i\^Ar, p{A) =N + k]. 
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Figure 3: A compact directed site animal (filled circles) on the triangular lattice with 
perimeter sites marked by open circles. 

The site case 

An animal is compact if tlie occupied sites in any given row are consecutive, i.e., tliere 
are no lioles in tlie animal (see Figure §). Obviously, removing interior sites from a 
compact animal can never reduce the perimeter. Therefore, the animals in An+i\An 
with minimal perimeter are compact. The minimal perimeter of a compact animal 
of height is + 2. This is proved by induction on A^. It is obviously true for 
A^ = 1 and one can easily see that by adding sites in row A^ + 1 to a compact animal 
of height A^ at least one more perimeter-site is added. We also note that there are 
at least two animals of height A^ with perimeter A^ + 2, namely a string of sites (one 
per row) running down either the left or right hand side of the lattice. This shows 
that A^ = A^ + 2. It is also clear that these two animals must be the ones that give 
rise to the first correction term dj^ Q = 2. What remains is to prove that there can be 
no more animals in An+i\An with perimeter A^ + 2. In order to do this we need a 
unique way of characterising the perimeter of compact animals of height A^. Introduce 
lines Rk {Lk) parallel to the right-hand (left-hand) egde starting from row k. Since the 
animal is compact all sites in A intersecting and are consecutive. The number 
of perimeter sites on the left-hand side of the animal is wi = max{k, (1 A ^ ^} 
because the last occupied site in line Lk has an unoccupied neighbour on L^. Similar 
arguments apply for the number of perimeter sites Wr on the right side of A. Finally 
we note that the only perimeter site not accounted for is the one lying vertically below 
the last site in Ljy and/or R^. So the perimeter is p{A) = wi + Wj. + 1. Furthermore 
if A G An+i\An then either wi or Wr (possibly both) has to equal A^. The animals 
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Figure 4: The two types of compact directed site animals with wi = 2 which contribute 
to the second correction term. 



with minimal perimeter + 2 are those with wi = 1 or Wr = 1, obviously there can 
be only two such animals, which completes our proof that dN,Q = 2. 

From Eq. ( |A.4| ) we get the second correction term 



where |.4.Ar,i| is the number of animals of height with a perimeter of length + 3. 
From the characterisation of compact animals derived above it follows that the animals 
in 1 are those with wi = 2 or Wr = 2. Obviously there is the same number of animals 
in each case so we can restrict ourselves to the case wi = 2, Wr = N. We are thus 
looking at animals restricted to the left-most two lines Li and L2 of the lattice and 
either LiHRn or L2f\RN has to be non-empty. The two types of animals are illustrated 
in Figure ^. If Li fl Rjq 7^ (Figure ^) then the first N sites of Li are occupied and 
1 < A; < consecutive sites along L2 are occupied; these k sites can be placed in 
N -k + l positions. If Li n i^jv = (Figure |b) and the first k sites (1 < A; < A^ - 1) 
of Li are occupied then the first j consecutive sites < j < A; of L2 may be empty. 
Combining these two contributions with those from = 2 we find 

/ N N-l \ 

\An,i\ = 2{Y,{N -k + l)+Y.{k + l)\ = 2N^ + 2N -2. 

\k=l k=l ) 

Since the number of sites in each of the two animals in A^^ is A^, Eq. ( [A.5| ) yields 
thus proving the empirical formula derived previously. 
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the third correction term. 



C2 

3 which contributes to 



Next we prove the formula for dN,2- From Eq. ( |A.4| ) we see that the third correction 
term is given by 



dN,2 = \A^,2\ - E (1^1 - 1) + E f '"^2 M • (A-6) 

In this case there are two distinctly different sets of animals in An,2, namely, compact 
animals with wi = 3 a.s pictured in Figure and animals formed from the compact 
animals of Figure ^ by removing consecutive sites from the second line of occupied 
sites leaving at least the first and last sites untouched. One easily sees that cutting 
such a 'hole' in these animals is the only way of increasing their perimeter by one site. 
From the animals in Figure |^ we get the following contributions 



TV k 

a: 2EE(^-'^ + l)('^-^ + l) 

k=l 1=1 

N-1 k N~l 

b: 2EEE(^-^-^ + i) 

k=l 1=1 m=l 

N-1 k 

Ci: 2EE(^ + 1) = 

k=l 1=1 
N-2 k N-k-1 
C2: 2EE E (^ + ^ + 1) = 

k=l 1=0 m=l 



—N"" + -N"^ + —N^ + -N. 



12 



12 



4 6 4 6 ' 



-N^ + - -N, 
3 3 ' 

6 3 6 3 



(A.7) 



The animals in Figure ^ account for animals with Li fl Rn ^ 0, those of 03 for 
animals with Li fl Rn = and L2 fl Rjy 7^ 0, and lastly those of ^ for animals where 
Li n Rn = and L2 fl Rn = 0- The contribution in each case is simply all the possible 
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configurations which leads to an animal of the specified kind. The sums in Eq. ([A.7|) 
should be self-evident. 

The animals in Figure El with a cut as described above yield the contributions 



N k-2 

a: 2EE(^-^ + l)(^-^-l)= I^^' - 6^' - I^^' + 6^' 

N-lk-2k-l-l 1111 

b: 2Y.T.T.ik-l-m)= T-N' - -N' - -N' + -N . (A.8) 

k=2 1=0 m=\ 1/ D iZ D 

In ) the piece in the second line has to have at least three sites (/c > 3) otherwise 

one could not cut out a hole of size / < k — 2. The k sites can be placed in (A^ — k + 1) 
positions and the hole can be cut in k — 2 — 1 + 1 = k — I — 1 places, which leads to the 
first sum. In case (b) there can be from 2 to — 1 sites in the first line (the sum over 
k) with an overlap ofO<m<k — 2 sites between the first line and the consecutive 
sites in the second line extending to the A^th row. Among the remaining k — m sites 
in the second line l<l<k — m — 1 are occupied and they can be placed mk — m — l 
positions, thus giving us the second sum. 



The second term in Eq. ( |A.6|) is the sum over | A| — 1 of the compact animals in Figure ^ 



and we find the two contributions: 



N A 1 

a: 2Y,{N -k + l){N + k-l) = -N^ + - -N, 

k=i 3 3 

7V-1 k 4 in 

b: 2j:T.iN + k-l-l)= -N'--^N + 2. (A.9) 

fc=l 1=0 

Finally the last term in Eq. ([A.6|) simply stems from the two animals in An,o and their 
contribution is 



^(^2^) ^ N^-3N + 2. (A.IO) 



By adding the contributions of Eqs. (|A.7| ), ( [A.8| ) and ( [A.lOp while subtracting those 
of Eq. (^) we get 



dN^2 = -N^ -N^ + -N"^ -2N + 2 = — (SAT^ - 12N^ + AN"^ - 24:N + 24) (A. 11) 
3 3 12 

in full agreement with the extrapolation formula listed in Table |I|, thus concluding the 
proof for dN,2- 
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Figure 6: A site-compact directed site-bond animal (filled circles and thick bonds) on 
the triangular lattice with possible perimeter sites marked by open circles. Some of the 
perimeter sites have only one possible incident bond (marked by double lines) and in those 
cases the bond can be present (the site is part of the perimeter) or absent (the edge is part 
of the perimeter). 

The site-bond case 

From the emperical extrapolation formulas it it clear that the site-bond case is very 
similar to the site case and only a few generalisations are necessary. Again we look at 
compact animals and the ones we shall call site- compact have the minimal perimeter. 
A site-compact animal is one in which, as before, all occupied sites and bonds in a row 
are consecutive and in addition all possible bonds to sites with more than one incident 
edge are present. Figure ^ shows such an animal. Clearly the perimeter of such an 
animal is equal to the perimeter of the identical site animal. Thus the animals with 
minimal perimeter have wi = 1 (or Wr = 1). Such animals consist of consecutive 
occupied sites down the left-hand side with most of the bonds emanating from these 
sites present. A few of the bonds can be either present or absent, namely, the bond 
from the top site pointing South-East and the bonds from the last site pointing South- 
West or South, though in this latter case at least one of the bonds has to be present. 
So all in all there are three possible bond configurations from the last site and two 
from the top site for a total of six possibilities. Taking into account the animals with 
Wr = I we have proved 

dN,o = 12 
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Figure 7: A compact directed bond animal (thick bonds) on the triangular lattice with 
perimeter bonds marked by open circles. 

The bond case 

The first correction term for tlie bond case, d^^o = 2Cn — 1, involve tfie Catalan 
numbers Cat which equal the first correction term for the square bond problem (Baxter 
and Guttmann 1988). Bousquet-Melou (1995) proved this result by noting that the 
square bond correction term arise from compact bond animals of directed height N. 
The first correction term for the triangular bond problem can be found by generalising 
the arguments from the square bond case. The first correction arise from compact 
animals constructed as follows. Choose two paths Ui and 002 consisting of bonds 
pointing only South and South- West starting from the origin and terminating at the 
same point on level N. The animal obtained by filling in all bonds between Ui and 
LJ2 has height N and perimeter 2N + 1. These animals are just the staircase animals 
which are enumerated by the Catalan numbers and give rise to the first square bond 
correction term. Obviously the set of animals bounded by paths consisting of South and 
South- East bonds also contribute to the first correction term. The animal consisting 
entirely of south bonds (a line of bonds down the center of the lattice) is the only 
animal included in both sets. The first correction term is exactly due to these 2Cn — 1 
'staircase animals' on the triangular lattice. 
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